Improved limits on short-wavelength gravitational waves from the cosmic microw^ave 

background 



o 

(N 



a\ 



Irene Sendra^'El and Tristan L. Smith^'^'0 

Fisika Teorikoa, Zientzia eta Teknologia Fakultatea, 

Euskal Herriko Unibertsitatea UPV/EHU, 644 Posta Kutxatila, 48080 Bilbao, Spain 

^Berkeley Center for Cosmological Physics, Department of Physics, 

University of California, Berkeley, CA, USA 94720 

* Institute for the Physics and Mathematics of the Universe (IPMU), University of Tokyo, Chiba 277-8582, Japan 

The cosmic microwave background (CMB) is affected by the total radiation density around the 
time of decoupling. At that epoch, neutrinos comprised a significant fraction of the radiative energy, 
but there could also be a contribution from primordial gravitational waves with frequencies greater 
than ~ 10~^^ Hz. If this cosmological gravitational wave background (CGWB) were produced 
under adiabatic initial conditions, its effects on the CMB and matter power spectrum would mimic 
massless non-interacting neutrinos. However, with homogenous initial conditions, as one might 
expect from certain models of inflation, pre big-bang models, phase transitions and other scenarios, 
the effect on the CMB would be distinct. We present updated observational bounds for both initial 
conditions using the latest CMB data at small scales from the South Pole Telescope (SPT) in 
combination with Wilkinson Microwave Anisotropy Probe (WMAP), current measurements of the 
baryon acoustic oscillations, and the Hubble parameter. With the inclusion of the data from SPT 
the adiabatic bound on the CGWB density is improved by a factor of 1.7 to 10®r2g„/i^ < 8.7 at 
the 95% confidence level (C. L.), with weak evidence in favor of an additional radiation component 
consistent with previous analyses. The constraint can be converted into an upper limit on the 
tension of horizon-sized cosmic strings that could generate this gravitational wave component, with 
G/x < 2 X 10~^ at 95% C. L., for string tension G/i. The homogeneous bound improves by a factor 
of 3.5 to lO^Jlgw/i^ < 1.0 at 95% C. L., with no evidence for such a component from current data. 
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I. INTRODUCTION 



A primordial gravitational wave background may have 
been generated by processes taking place in the early Uni- 
verse. Given the extremely small cross section of gravita- 
tional waves (GWs), they can probe deep into the early 
Universe, and provide a unique window to explore its evo- 
lution. The exact mechanisms of production of primor- 
dial GWs are still under investigation. There are many 
theoretical models that predict them: not only through 
quantum fluctuations during inflation ^E] , but also from 
cosmic strings [IHS] , causal mechanisms from phase tran- 
sitions PUHH] , ekpyrotic models [121 [H] or pre big-bang 
theories [T71 [TH]. Such backgrounds will involve GWs 
with wavelengths that extend up to our present cosmo- 
logical horizon, giving a lowest observable frequency limit 
of- 10-1^- 10-16 Hz. 

A measure of the amplitude of the cosmological 
gravitational-wave background (CGWB) at low frequen- 
cies can be obtained from constraining a possible tensor- 
mode contribution to the large-scale temperature and 
polarization fluctuations in the cosmic microwave back- 
ground (CMB) [T9H22] . Recent results from the Wilkin- 
son Microwave Anisotropy Probe (WMAP) satellite limit 
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the amplitude of tensor fluctuations, quantified by the 
tensor-to-scalar ratio, to r < 0.20 at the 95% confldcnce 
level (C. L.) [23], which translates to Vl^^h? < lO^^^ at 
frequencies ^ 10" ^^ — 10" ^^ Hz. At higher frequencies, 
larger than — 10^^° Hz, a bound on the CGWB can be 
obtained through big-bang nucleosynthesis (BBN). Grav- 
itational waves at these frequencies would contribute to 
the total radiative energy density at the time of nucle- 
osynthesis, thus a measure of the light-element abun- 
dances can place a constraint on the CGWB. Current 
observations set an upper limit on the CGWB energy 
density at this frequency range of lO^fig^/i^ < 8.0 
[24j . At the same time, pulsars act as natural grav- 
itational wave detectors |25j . strongly constraining the 
amplitude of the GW with the narrow range of frequen- 
cies, / ~ 10-9 - 10-8 Hz, to 10Sf]g«/i2 ^ ;lq-2 gg 

Large scale interferometers for gravitational wave detec- 
tion (e.g., LIGO [27], LISA ^, VIRGO [29]) are also 
looking for gravitational wave signals. A recent bound 
of lO^rJgw^^ < 6.9 has been obtained at 10^ Hz from 
the Laser Interferometer Gravitational Wave Observa- 
tory (LIGO) [30]. 

In this work we estimate constraints on the CGWB at 
frequencies larger than — 10- ^^ Hz through their effect 
on the angular power spectrum of the CMB. This compo- 
nent behaves as a non-interacting relativistic fluid, and 
so would modify the CMB power spectrum in the similar 
way as adding extra neutrino species [31] . Given the — 2ct 
hint of a neutrino number excess seen by both the South 



Pole Telescope (SPT) and Atacama Cosmology Telescope 
(ACT) experiments [32l [33], it is timely to consider al- 
ternatives beyond the standard three neutrino species. 
Constraints on this GW background were first presented 
in [31,, and updated with WMAP 7-year data in [25] , 
We now make use of the latest CMB data from WMAP 
and SPT, which currently give the tightest constraints on 
the number of effective neutrinos [35] . We then use the 
improved limits to bound the tension of cosmic strings 
that could generate this gravitational wave background. 
As in [3T] , we do not restrict our study to the case where 
the GW where produced under adiabatic initial condi- 
tions, but also study homogeneous initial conditions. In 
Section \n\ we introduce the theoretical background, in 
Section |III| we present the observational data used, and 
finally, in Sections IV and[V]we show the results obtained 
and discuss our conclusions. 



II. THEORETICAL BACKGROUND 

The effective number of neutrino species, N^ff, repre- 
sents the energy density stored in relativistic components 
(radiation) as 



Prad ^ P-t + Pu + Px 
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with p-y, p^, and px the energy densities of photons, neu- 
trinos, and possible extra radiation components, respec- 
tively. In this expression, all non-photon energy density 
is expressed in terms of an effective number of neutrino 
species, N^s- Any additional radiation density, px, would 
then typically correspond to a non-integer value of N^s. 
The largest effect of increasing the radiation density on 
the CMB comes from a decrease in redshift of matter- 
radiation equality, Zeq, which follows from the relation 



N^s = 3.04 + 7.44 



VLmh^ 3139 

0.1308 1 + Zeq 



1 



(2) 



for matter density ilmh"^ [HI]. An increase in N^a also 
affects the acoustic oscillations in the primordial photon- 
baryon plasma, leading to additional damping and a 
phase shift in the acoustic peak positions of the CMB (see 
[55]). As noted in [3S], the matter density can be better 
constrained by combining observations of the CMB with 
late-time distance measurements from baryon acoustic 
oscillation (BAG) data, and a measurement of the Hub- 
ble parameter, Hq. These allow the equality redshift to 
be better measured, improving constraints from small- 
scale CMB data on iVeff. 

If we assume that there are three standard neutrino 
flavors, A^cff = 3.046 (with the small correction due to 
finite temperature QED effects and neutrino flavor mix- 
ing |36j). and any measured excess would imply an extra 
relativistic component. Current small-scale CMB data 
show a slight preference for Ncs > 3 at 95% C. L., 



with A^off = 4.56 ± 1.5 measured from the Atacama Cos- 
mology Telescope and A^eff = 3.86 ± 0.84 with South 
Pole Telescope data [35] in combination with WMAP 
data, measurements of the BAG, and Ho [33;. Combin- 
ing both ACT and SPT (along with WMAP and other 
probes of large-scale structure) increases the significance 
to A^cff = 4.0 ± 0.58 [37]. As proposed in f3T], such an 
excess could be interpreted as a contribution of gravita- 
tional waves to the radiation density, instead of neutri- 
nos; the effects on the CMB and the matter power spectra 
are similar in both cases. Since gravitational waves with 
wavelengths shorter than the sound horizon at decou- 
pling behave as free-streaming massless particles, a limit 
on additional relativistic radiation can be translated into 
an upper limit on the energy density of the CGWB for 
frequencies larger than ~ 10""'^^ Hz. Gne can relate the 
effective 'number of gravitational-wave' degrees of free- 



dom, Afg 
!35': 



A^ni 



3.046, to the CGWB energy density 
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d(ln/)/i2^3^(/) = 5.6 X lO-e^Vgw. (3) 



The different processes which may produce a CGWB 
lead to different initial conditions for perturbations in 
the CGWB fluid. Adiabatic initial conditions may arise 
from the incoherent superposition of cusp bursts from 
a network of cosmic strings or superstrings (see [T] [8] 
for further details). Bounds on the energy density of 
this radiation can then be translated into an upper limit 
on the string tension, G^, for a given model [IHEl 135] . 
Following the procedure given in [8] for the analytical 
approximation of the CGWB produced by strings, we 
can relate our limit on the GW energy density to a bound 
on Gp and the string reconnection probability, p, which 
is 1 for cosmic strings and < 1 for superstrings. String 
loops which are smaller than the causal horizon produce 
a CGWB with 



^gw(/) ~ 5 X IQ-'^Gp/p 



(4) 



whereas the CGWB produced by horizon sized loops 
scale as 



^gw(/) ~ 3.2 X IQ-^^/Gp/p. 



(5) 



An alternative scenario, noted in |31j . may arise when 
the CGWB is produced by quantum fluctuations during a 
period of inflation. In this case, the perturbations may be 
non-adiabatic, differing from those of other species pro- 
duced though the decay of the inflation. We re-examine 
the dependence on initial conditions by also considering 
'homogeneous' initial conditions, as in [3T]. These have 
no initial primordial perturbations in the CGWB energy 
density in the Newtonian gauge, thus the curvature per- 
turbation disappears when the CGWB energy density 
dominates, and, in addition, it approaches the adiabatic 
case when the CGWB energy density vanishes. 



III. DATA AND METHODOLOGY 

We use the latest measurements of the CMB power 
spectrum from the WMAP 7-year data release [40] and 
SPT f32] , which extend the previous set of data to higher 
multipole moments (^max ~ 3000). These data consider- 
ably improve constraints on iVoff , which will be translated 
into tighter bounds on flg^h^. To illustrate the sensitiv- 
ity of these new data to changes in iVoff , in Fig. [T] we 
show the difference between a model with N„^ = and 



N^- 



1 with the other cosmological parameters chosen 



so that the two models fit the WMAP-7 data equally 
well. It is clear that the data from SPT covers the range 
where the power spectrum is sensitive to a variation of 
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We estimate the bounds on iVgw by varying the number Figure 1: The difference between the temperature power spec- 



of effective gravitational-wave degrees of freedom, iVg. 
imposing a prior that N^^ > and assuming the stan- 
dard three non-interacting massless neutrino species and 
by marginalizing over the standard six ACDM parame- 
ters: the baryon energy density in units of the critical 
energy density, fl^h^, the energy density in cold dark 
matter in units of the critical energy density, fluuh'^, 
the angular size of the first acoustic peak 9, the opti- 
cal depth to the surface of last scattering, r, the scalar 
spectral-index, Us , and the amplitude of the scalar power- 
spectrum, As ■ To include the SPT data we follow the pre- 
scription presented in |32j , marginalizing over three addi- 
tional parameters describing power from foregrounds and 
the Sunyaev-Zel'dovich effect. We use a Markov Chain 
Monte Carlo (MCMC) method to estimate the probabil- 
ity distribution, using the CosmoMC^ software package, 
see [12135] with a modified version of CAMB^ |43j which 
includes the homogeneous initial conditions for CGWB 
perturbations in the conformal Newtonian gauge. We 
combine these CMB measurements with distance mea- 
surements from Baryon Acoustic Oscillation (BAO) data 
measured at z = 0.2 and z = 0.35 from SDSS and 
2dFGRS presented in [l^. We also include a prior of 
Ho = 74.2 ±3.6 km/(s Mpc) on the Hubble constant [45] . 
The CMB temperature and polarization spectra are ob- 
tained with CAMB for each cosmology, and then passed 
to CosmoMC to obtain the joint likehhood of WMAP, 
SPT, BAO, and Hq, as described in &. 



with the other cosmological parameters chosen so that the 
two models fit the WMAP-7 data equally well. The thick-line 
model has adiabatic initial conditions, the dot-dashed homo- 
geneous initial conditions. The angular range of the SPT and 
WMAP 7-year CMB data are shown. 
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Figure 2: Marginalized likelihoods for A'^gw and CGWB en- 
ergy density, Qg„h^{xlO~^), for adiabatic (left) and homoge- 
neous (right) initial conditions, where the horizontal line and 
the colored/shaded region indicate the 95 % C.L. limit. 



IV. RESULTS 



As shown in Fig. [2] and Table [T] the new combination 
of CMB data gives an upper limit of iVgw < 1.56, or 
lO'^rigw^^ < 8.71, at 95% C. L. for adiabatic initial con- 
ditions which improves the previous upper bounds on 



^ http://cosmologist.info/cosmomc/ 

^ http://camb.info/ 

^ Likelihood code is available at http://lambda.gsfc.nasa.gov/ 



Vtg-vjh? given in [23j by a factor of 1.7. A non-zero rel- 
ativistic excess is weakly favoured at 95% C. L., with 
the best-fit value A^g^ = 0.82 ± 0.36. Homogeneous ini- 
tial conditions for GWs are now more strongly excluded, 
with Agw < 0.18, lO^Jlgwft.^ < 1.02, an improvement on 
the previous constraint in [9i by a factor of 3.5. In this 
case the degeneracy between the effects of CGWB and 
neutrinos is broken, giving us the stronger constraint. 
These bounds are now competitive with those given by 
BBN, 10^f7gw/i^ < 8.0 ^24J for / > 10"^^" Hz and LIGO, 
lO^f^gw/i^ < 6.9 at / w'lO^ Hz [5^, but with the benefit 
that this upper bound can be extended to frequencies as 



Table I: Upper limits on A'^gw and ilgwft^ at 95% C.L for the 
adiabatic and homogeneous primordial initial conditions. 



Adiabatic 95% upper limits 

Adiabatic 68% range 

Homogeneous 95% upper limits 



JVgw Wilg^h^ 



< 1.56 < 8.71 
0.82 ±0.36 4.8 ±0.2 

< 0.18 < 1.02 



low as IQ-i^ Hz. 

Our constraint to the integrated CGWB energy den- 
sity with adiabatic initial conditions can be translated 
into a constraint on cosmic string networks, as shown 
in Fig. [3] In order to do this we use Eq. (pi) and inte- 
grate the spectra, given in Eqs. Q and (Is]). The lower- 
bound to the CGWB spectrum produced by the string 
network is fmin ^ 3.6 x 10^^^/(G/i) Hz for horizon-sized 
string loops and fmin ^ cn^^Zcq Hq for sub-horizon sized 
string loops where Zoq — 3400 is the redshift at matter- 
radiation equality ^. The upper-bound is given by the 
horizon size at the time of the phase-transition which pro- 
duced the network, f^^^ ~ a-^H{Tpt) = a-^T^jMpi = 
a~^MpiGfj., where Tpt = y/GjxMpx is the temperature of 
the phase transition, Afpi — I/^/SttG is the Planck mass, 
and a = 1 for horizon-sized string loops ^. With this 
we find that for subhorizon-sized loops our constraint be- 
comes, 






-^cq 



and for horizon-sized loops we have, 



In 



(GM)'Mpi 
3.6 X 10-18 Hz 



< 5 X 10"^ 



(6) 



(7) 



For p = 1 we obtain Gfi < 2 x 10^^ for horizon-sized 
string loops and G^jl < 2.5 x 10^^ for subhorizon loops. 
String tensions this large have already been excluded 
through directly limiting a string contribution to the 
CMB power spectrum [33l|46H49], Gfi < 10"*^ - 10"^. 
However, for lower reconnection probability, e.g. p — 0.1, 
the A'gw bound tightens to G^ < 2.4 x 10~^ for horizon- 
sized string loops. Further study is required to investi- 
gate constraints considering both the string contribution 
to the CMB power spectra, and additional relativistic 
degrees of freedom. 



V. CONCLUSIONS 

We have used recent CMB data at small angular scales 
from SPT, in combination with data from WMAP-7 year 
and the latest measures of Hq and BAO, to constrain 
a possible cosmological gravitational wave background 
with frequencies greater than 10~i^ Hz. Recent mea- 
surements of the Silk damping tail from SPT improves 
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Figure 3: The shaded region shows the 95% C. L. region 
which is allowed by constraints on the CGWB from the CMB 
in the p — G/x plane for horizon-sized (solid) , and subhorizon- 
sized (dashed) cosmic string and superstring models with 
string tension G\i and reconnection probability p, given a GW 
energy density level figw/i^ = 8.7 x 10"''. For p = 1, the CMB 



power spectrum limits G/U < 10 



10" 



limits on the CGWB density by about a factor of two. We 
note that the inclusion of additional measurements of the 
small-scale CMB, such as from the Atacama Cosmology 
Telescope, do not significantly improve the upper-limits 
presented here j57] . 

In agreement with other analyses of similar data we 
find that a non-zero relativistic energy density with adi- 
abatic initial conditions is preferred at the 95% level 
[32l [33l l37] . This additional energy density could be at- 
tributed to gravitational waves with frequencies greater 
than ~ 10^^ Hz, or some other relativistic species (see, 
e.g., [SZ])- If we attribute this extra relativistic en- 
ergy density to gravitational waves produced by a cos- 
mic string network it is interesting to note that for p = 1 
the upper limits to the CGWB presented here are com- 
petitive with direct CMB constraints on the the string 
tension, both now requiring G[i < 10"^ — 10"^. 

We have also explored how constraints on the CGWB 
depend on the initial conditions. Besides the standard 
adiabatic initial conditions, we placed constraints on the 
CGWB energy with 'homogenous' initial conditions ^. 
The additional information on the small-scale CMB from 
SPT improves upon the previous upper-limit by more 
than a factor of 3.5. As opposed to an additional adi- 
abatic relativistic component, the data does not prefer 
any additional homogeneous relativistic energy density. 

Finally, the situation will be greatly improved with 
near-future CMB observations. For example, the Planck 
satellite will increase the sensitivity of these constraints 
by a factor of 5 for adiabatic initial conditions and a 
factor of 2 for homogeneous initial conditions ^ . 
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